Low temperature transport in the XXZ model 
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We present evidence suggesting that spin transport in the gap less phase of the S=l/2 XXZ 
model is ballistic rather than diffusive. We map the model onto a spinless fermion model whose 
charge stiffness determines the spin transport of the original model. By means of exact numerical 
diagonalisation and finite size scaling we study both the stiffness and the level statistics. We show 
that the stiffness is non-zero at low temperatures so that the transport is ballistic. Our results 
suggest that the non-zero stiffness is due to the fact that even in the presence of Umklapp scattering 
a non-zero fraction of states remain degenerate in the thermodynamic limit. 
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The problem of transport in a non-disordered inter- 
acting many particle system is one of the oldest unsolved 
problems in solid state physics. A particular case which 
has attracted recent attention is spin diffusion in one di- 
mensional spin systems. Recent experiments indicated 
that S=l chains with the gap in the excitation spectrum 
display diffusive behavior in reasonable agreement 
with theoretical work Q relating diffusion to classical 
scattering of excitations near the gap edge. On the other 
hand, measurements on gapless S=l/2 chains ||] show a 
different behavior. The authors fit their data with a dif- 
fusion constant which is much larger than either the value 
found experimentally for the S=l chains or the value 
Vs ~ Jy/2TTS{S + l)/3 expected from classical consid- 
erations [Q. We believe that the measured value for Vs 
in the S=l/2 system is so large that it implies that the 
diffusion is not an intrinsic property of an S=l/2 spin 
system but is due to a weak coupling to other degrees of 
freedom (for example to phonons |^). 

A theoretical analysis based on a continuum limit Lut- 
tinger liquid representation H] suggested that the dif- 
fusion constant was associated with Umklapp operators 
and was finite but exponentially large in the Umklapp 
gap. Damle and Sachdev |^ provided a detailed analysis 
indicating diffusive behavior for the S=l chain, however, 
their work suggests that in general low energy excitations 
of gapped systems display diffusive behavior even if the 
effective low energy theory describing these excitations is 
integrable. On the contrary, Zotos and collaborators have 
argued that integrable models exhibit ballistic transport 
while non-integrable models are diffusive Recently 
Fabricius and McCoy have shown that numerical com- 
putations of the long time behavior of the correlation 
functions of the S=l/2 XXZ chain in the T — oo limit 
are consistent with ballistic transport if the model has a 
XY anisotropy but suggest that at the isotropic Heisen- 
berg point the transport is not ballistic. Very recently, 
Monte Carlo and Bethe-ansatz analyses of the stiff- 
ness of related models have appeared, reporting similar 
conclusions. 



In this paper we approach the question in a different 
way. We use exact numerical diagonalisation and finite 
size scaling to study the spin stiffness of the XXZ model 
and the behavior of the energy levels which leads to a 
non-zero stiffness. 

The XXZ model is defined by the Hamiltonian 
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At T = the gapless phase — 1 < A < 1 of this model is 
characterized by a non-zero spin stiffness Ds [pd| (in the 
gapped phase |A| > 1 the stiffness is zero). The ques- 
tion of interest here is whether at T > the stiffness 
Vs remains non-zero, implying ballistic transport, or Vs 
vanishes, implying non-ballistic (and perhaps diffusive) 
transport. 

The XXZ model is equivalent via the Jordan- Wigner 
transformation to the spinless fermion model, 
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{clcn+i + h.C.) 
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In this mapping the fermion density-density correlation 
function represents the — correlator and therefore 
the description of fermion transport directly translates 
into the spin language. In particular the real part of the 
frequency dependent conductivity may be written as 



TZe a{uj) = 2t:Dc5{uj) + Uregi^) 
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defining the charge stiffness Dc. If I?c 7^ the model has 
infinite conductivity whereas if Dc = one has either a 
normal conductor {Dc — 0, Urcgi'.^J ^ 0) > 0) with diffu- 
sive transport or an ideal insulator [Dc = 0, (Jrcg{^ ^ 
0) = 0). In the spin language Dc becomes the spin 
stiffness, which, if non-zero, corresponds to ballistic spin 
transport, whereas if it is equal to zero then the long 
time relaxation at finite temperatures is expected to be 
diffusive. 
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As was first noted by Kohn |T^, in systems with pe- 
riodic boundary conditions the stiffness at T = can be 
related to the response of the ground state energy Eq to 
a magnetic flux </>, which modifies the hopping term in 
the Hamiltonian Eq. by the usual Peierls phase factor 
t ^ iexp(±z(/)/L), so that = {L/2)d'^Eo/dcj)'^[cj) ^ 0), 
where L is the system size. This phase factor can be 
absorbed into twisted boundary conditions for the wave 
functions. 

Kohn's method has been recently generalized to finite 
temperatures |7j, 
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where Z is the partition function of the system and n 
labels exact eigenstates. 

We further rewrite Eq. fl) as Dc = Di + D2, with 
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calculates the eigenvalues through a series of orthogonal 
transformations with accuracy given by machine preci- 
sion. Our choices of cj) lead to ^ 10"^ accuracy for the 
derivative values. The size of matrices that could be diag- 
onalised by the routine is limited by computer memory; 
for a X matrix it requires ~ SA"^ bytes of storage 
space. With the availiable computer memory of about 
360 MB we can diagonalise matrices up to A^ = 7000 
which limited us to chain sizes L < 14. 




The advantage of this representation is that it sepa- 
rates Dc into a thermodynamic part (depending only 
on derivatives of Z) and a positive part, depending on 
current-carrying (j„„ dEn/d4>) states. The thermody- 
namic part may be seen to give no contribution to the 
charge stiffness at T > for example at low tem- 
peratures T <C J 



Di{T > 2tt /vRL) = iTcxp(- 
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with i?2 = (1 - (l/7r)cos"i A)/27r. Thus at T > 
any non-zero Dc must be due to D2 which essentially 
counts the number of thermally accessible current carry- 
ing states. 

Time reversal invariance implies that for a non- 
degenerate state dEnld(j){<t> — > 0) = 0. Current carrying 
states occur in degenerate pairs which which are split by 
the application of magnetic flux. A sufficient condition 
for a non-vanishing Dc is to have a non-zero fraction of 
current carrying states with dEn/d(f> ^ In what 

follows, we investigate Dc and the statistics of the current 
carrying states numerically. 

To compute the eigenvalues we notice that for any fi- 
nite size chain the Heisenberg Hamiltonian Eq. is just 
a hcrmitian matrix. Our strategy is to construct this ma- 
trix for (/) = 1 X 10-^,2 X 10-"*, 3 x 10-'*,4 x lO""* and 
use the exact diagonalisation to obtain the eigenvalues 
and compute the derivatives. We use the standard QL 
routine from the Numerical Recipes package which 



Fig. 1. D2{T) for different system sizes for A = 0.4. 

The result of the calculation is presented in Fig. 1. 
For all system sizes we found D2 to be non-zero. At 
small temperatures the value of D2 appears to grow with 
system size (especially the peak value). At large tem- 
peratures all eigenvalues become involved in the sum Eq. 
(^ and the temperature dependence is defined by the 
prefactor 1/T, while the value of D2 decreases with the 
system size. 

To investigate the finite size scaling in Fig. 2 we plot 
{D2T)oo — lim D2T versus the inverse system size for 

T — >oo 

different values of the interaction. The symbols repre- 
sent the actual data points and the best fit lines are 
continued to the infinite size (1/i = 0). We are un- 
aware of theoretical results for the large L behavior of 
{D2T)oo] our numerical results are consistent with the 
ansatz {D2T)oo{L) = A + B/L + ... with A,B depending 
on the interaction, but with A always positive for A < 1. 
For small A the best fit line is flat and the fit using four 
largest sizes is excellent (error estimate for parameter A 
is 0.7%). For A = 0.6 we find A « 0.076, B w 0.32 with 
3.6% error. At the isotropic point A = 1 the best straight 
line fit yields A w 0.029, B sa 0.46 but with rather larger 
11% error leading us to question whether we have as- 
sumed the correct functional form. We note, however, 
that fits to the form {D2T)oo{L) — C/L^ lead to even 
larger errors, so the hypothesis {D2T)(L — > 00) — > is 
inconsistent with our data. 
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Fig. 2. [D2T) plotted against inverse system size at 
T = 50. 

The scaling of D2 can be expressed in terms of the 
size dependence of the current carried by a typical ex- 
cited state. For the free case a typical state contains a 
total number P ~ L of fermions excited above both left 
and right Fermi points and has a \/P imbalance between 
the left and right movers, producing a non-zero current 

dE^/dci> ~ i/VI. 



mensurate with a reciprocal lattice vector, they cannot be 
mixed by the Umklapp interaction. The interaction can 
mix a given current-carrying state with another current- 
carrying states changing the value of the total current 
carried. To analyse these effects we plot in Fig. 3 the 
fraction of states with dEn/d4>{4> 0) = as a function 
of system size for the XXZ model with varying interac- 
tion strength A (sohd symbols). One sees that adding an 
interaction sharply increases the fraction of non-current- 
carrying states, but this fraction remains small and de- 
creases with system size. 
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Fig 3. Fraction of states with zero current 

The interaction affects current carrying states in two 
ways. As we noted above these states come in degenerate 
pairs. If the momentum of two degenerate states differs 
exactly by a reciprocal lattice vector, the these states 
will be mixed by the Umklapp interaction term which 
will destroy the degeneracy. Consequently these states 
will no longer carry current. On the other hand if two 
degenerate states differ by a momentum which is incom- 



Fig 4. Histogram of current values; lines for different 
A are indistinguishable 

We now consider the statistical distribution of the cur- 
rents carried by the eigenstates. We show in Fig. 4 a 
histogram of {dEn/d(j))'^ values for L = 14 and all previ- 
ously considered values of A. For each A the a;- axis has 
been scaled so a; = 10 corresponds to {dEn/dcj))'^ equal 
to the average value. The data have been grouped in to 
bins of width 0.1 of the average {dEn/d(j>)'^ and the y- 
axis has been scaled so that y = 1 represents the number 
of states in bin l.With this choice of scaling the distribu- 
tions for different A are indistinguishable: the interaction 
does not change the shape of the distribution, but merely 
reduces the avarage value of {dEn/ dcff' . 

We now consider the effect of an interaction that spoils 
the integrability of the XXZ model by adding the next- 
nearest neighbor interaction V ^ SlSf^2 to the Hamilto- 

i 

nian Eq. (yj). As shown in Fig. 3 this term lifts more de- 
generacies, so that more states carry zero current. How- 
ever the relative change in the fraction of these states is 
small and the size dependence is similar to the integrable 
case. 

The effect of the non-integrable interaction is mostly 
to reduce the values of the current carried by the remain- 
ing degenerate states as is illustrated on Fig. 5. At least 
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for large V, {D2T){L oo) 0, implying that the cur- 
rent carried by a typical state, although not 0, is much 
less than 1/VL and is presumably o(l/L). Interesting, 
the effect seems not to occur for F < A. For A = 
Fig. 5 shows clearly that even V ~ 0.02 leads to D2T 
which vanishes as L 00, whereas for A = 0.2 large 
V = 0.62 leads to a vanishing {D2T){L 00), while the 
effect oi V = 0.02 is much smaller than for A = and 
our data are consistent with a non-zero {D2T){L 00). 
Our system sizes are too small to allow us to make a def- 
inite statement about {D2T){L 00) for V < A, but 
clearly the relative size of the effect of the non-integrable 
interaction depends strongly on the ratio V/A. 
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the XXZ model. However, the degeneracy of the XXZ 
model is due to its integrability and it was earlier argued 
that the anomalous transport is due to the large num- 
ber of conservation laws which characterise an integrable 
model. We believe that among all those conservation 
laws there should be only one that is responsible for the 
degeneracies and therefore for the anomalous transport. 
Based on our results for the non-integrable model with 
the next-nearest neighbor interaction we infer that this 
conservation law is probably not specific to integrable 
models. The next-nearest neighbor interaction destroys 
the integrability but still leaves a macroscopic fraction 
of degenerate (current-carrying) states. In the integrable 
models (except possibly for A = 1) and possibly in the 
non-integrable models (for y < A) the typical current 
carried by these states is ^ 1/VL leading to non-zero 
stiffness as L ^ 00. For V > A we see the typical cur- 
rent is much smaller, leading to vanishing stiffness and 
speculate that there is a length scale f (V, A) such that 
for system sizes L < ^ the current carried by a typical 
state scales as 1/ \/L but that for L > ^ the current scales 

as r^/yi. 
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in the problem and for helpful discussions and to Ian Af- 
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paper. A. J. M thanks the Physics Department of Rut- 
gers University for hospitahty and NSF DMR 9707701 
for support. 
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Fig 5. D2T{1/L) comparing integrable and non- 
integrable cases; the lines are guides to the eye. 

In conclusion we have studied the spin transport in the 
Heisenberg model by calculating the finite temperature 
stiffness for small system sizes. The data presented in 
Fig. 2 show that for the available sizes D2 is greater than 
zero and extrapolates to a non-zero value in the ther- 
modynamic limit. At small A this extrapolation seems 
unambiguous in agreement with the perturbation theory 
results which are valid for small A where the Umk- 
lapp is irrelevant in the renormalisation group sence. For 
larger A the data seems to slowly decrease with size ex- 
trapolating to some small but non-zero value as 1/L —> 0. 
For the isotropic Heisenberg point (A = 1) the considered 
sizes are too small to make conclusive predictions about 
the thermodynamic limit behavior perhaps because the 
Umklapp becomes marginal at this point. 

The non-zero stiffness at T > means that the model 
posess a macroscopic fraction of current carrying states. 
These states come in degenerate pairs, so the stiffness 
essentially is a measure of non-trivial degeneracy of the 
energy eigenvalues. We argue that this degeneracy is the 
underlying physical reason for the ballistic transport in 
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